Could strange stars be in the color-flavor-locked phase: Tested by their thermal evolutions 
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The thermal evolution of strange stars in both normal and color-flavor-locked (CFL) phases are investigated 
together with the evolutions of the stellar rotation and the r-mode instability. The heating effects due to the 
deconfinement transition of the stellar crust and the dissipation of the r-modes are considered. As a result, the 
cooling of the stars in the normal phase is found to be not very different from the standard one. In contrast, for 
the stars in the CFL phase, a big bump during the first hundred years and a steep decay (~7% in ten years) at the 
ages of ~ 10 4 ~ 6 yrs are predicted in their thermal evolution curves. These unique features provide an effective 
observational test for determining whether or not the CFL phase is reached in strange stars. This thermal test 
method is independent of and complementary to the rotational test method, which is a direct consequence of the 
r-mode instability [see J. Madsen, Phys. Rev. Lett. 85, 10 (2000)]. 

PACS numbers: 97.60Jd, 04.40.Dg, 21.65.Qr, 97.60.Gb 



I. INTRODUCTION 

Strange quark matter (SQM), made up of roughly equal 
numbers of up, down, and strange quarks, could be the abso- 
lute ground state of strong interaction [1-3], the most impor- 
tant astrophysical consequence of which is that compact stars 
could have a quark matter core (i.e., hybrid star) or even con- 
sist purely of SQM (i.e., strange star, SS; [4, 5]). Furthermore, 
phenomenological and microscopic studies suggest that SQM 
at a sufficiently high density could undergo a phase transition 
into a color superconducting state, the typical cases of which 
are the two-flavor color superconductivity (2SC) and color- 
flavor-locked (CFL) phases [6, 7]. In the following content, 
we adopt the abbreviations 2SSs and CSSs for the two types 
of SSs in the two different color superconducting phases, re- 
spectively, and NSSs for normal SSs. 

Madsen [8] suggested that the SS models can be con- 
strained by the most rapidly rotating pulsars in low mass X- 
ray binaries (LMXBs), because the rotation of compact stars 
should be suppressed by r-mode instability, which arises due 
to the action of the Coriolis force with positive feedback in- 
creasing gravitational radiation (GR) [9? ]. As a result, it 
is found that the CSS model could be not permitted by the 
LMXB data because of the too weak viscosities of CSSs, 
which cannot inhibit the r-modes effectively. In this paper 
we suggest something different from the rotational test to the 
models — a thermal test method, i.e, probing the interior of 
SSs by the temperature observations of compact stars. The 
most essential point of our work is to carefully investigate the 
influence of the r-mode instability on the thermal evolution of 
SSs. Therefore, our independent thermal test can be comple- 
mentary to the direct rotational test. 

The cooling behaviors of SSs have already been extensively 
investigated, ever since the emergence of the SS model. In 
the earliest works, it was accepted that the surface temper- 
ature of SSs should be significantly lower than traditional 
neutron stars at the same age due to the quark direct Urea 



(QDU) processes llTl1 - tl4ll . However, by considering that the 
electron fraction of SQM could be small or even vanish and 
thus the QDU processes could be switched off, Schaab et al. 
11151 [lal proposed that the cooling of SSs dominated by the 
quark-modified Urea and quark bremsstrahlung processes can 
be slower than neutron stars with standard cooling. The pos- 
sible existence of the color superconductivity makes the prob- 
lem more complicated. It has been shown that the cooling 
of 2SSs is compatible with existing x-ray data but that CSSs 
cool down too rapidly, which completely disagrees with the 
observations [17]. 1 

Besides the stellar core consisting of SQM, SSs could also 
sustain a tiny nuclear crust with a maximum density below 
neutron drip (~ 4 x 10 11 g ctrT 3 ) 0, [H, [H] . The mass of the 
crust as a function of the spin frequency v can be expressed as 
fSJl 



M c = M° c (l + 0.24 v\ + 0.1 6vij) , 



(1) 



where v 3 = v/10 3 Hz, and M® ~ 10" 5 M o is the mass of the 
crust in the static case. With the spin down of a SS, the star 
shrinks and the density at the bottom of the crust exceeds neu- 
tron drip. As a result, the most inner part of the crust falls 
into the quark core and dissolves from baryons into quarks. 
During this process the released binding energy can heat the 
star. Therefore, Yuan and Zhang [21] proposed a so-called 
deconfinement (DC) heating mechanism, which can delay the 
cooling of SSs significantly. Especially for CSSs, the stellar 
thermal evolution could be completely determined by the DC 
heating B22I1 . As analyzed, the temporal dependence of the 
heating power should be determined by the spin-down history. 
In B21I1 and B22I1 . as usual, the spin down of SSs is consid- 
ered to be due only to the magnetic dipole radiation (MDR). 
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1 In the CFL case, both the neutrino emission and the quark thermal capac- 
ity are blocked off, whereas the cooling due to the photon emission from 
the stellar surface is still effective. Therefore, an extremely large temporal 
derivative of temperature can be obtained by dividing the photon luminos- 
ity by the tiny heat capacity, which leads to the rapid cooling. See the grey 
band in the right panel of Fig. [5] 



However, as aforementioned, the spin of young compact stars 
could also be limited by r-mode instabilities. In other words, 
if the r-mode amplitude is large enough, the early spin down 
of SSs could be dominated by GR rather than MDR [H . This 
could significantly influence the process of the release of the 
crust-binding energy. On the other hand, as an extra heating 
effect, the shear viscous dissipation of r-modes can also de- 
posit energy into the thermal state of the stars [24-26]. The 
main purpose of this paper is to test the CSS model by in- 
vestigating its thermal evolutions with the above-mentioned 
heating effects. 

In the next section we briefly introduce the thermal evolu- 
tion equation of SSs, where both DC heating and the heating 
due to viscous dissipation of r-modes (RM heating hereafter) 
are involved. The spin evolution of SSs involving both MDR 
and GR brakings is described in Sec. [HI] which is essential 
to solve the thermal evolution equation. The calculated re- 
sults and some discussions are provided in Sees. |IV|and[Vl 
respectively. 



II. THERMAL EVOLUTION EQUATION 

For simplicity, in this paper, the temperature in the inte- 
rior (the core) of a SS is assumed to distribute nearly uni- 
formly, and the gradient from the interior temperature to the 
surface temperature is considered to be mainly caused by the 
thin crust. Following 02711 . the relationship between the sur- 
face temperature T s and the interior temperature T can be ex- 
pressed by T s = 3.08 x 10 6 g' /4 4 r°- 55 , where T g = T/10 9 K 
and g s 14 is the proper surface gravity of the star in units of 
10 14 cm s~ 2 . As a result, the luminosity of the surface pho- 
ton emission of the star can be calculated to L y - AnR 2 crT 4 s , 
where R is the stellar radius and <x is the Stefan-Boltzmann 
constant. Of course, as usually considered, neutrino emission 
could play a much more important role in the cooling of a 
young NSS. 

Denoting the neutrino luminosity by L v and the powers of 
the DC and RM heatings by Hoc and //rm, the thermal evo- 
lution of a SS can be solved from 



C- 



JT 

Hi 



L v + Hdh + H\ 



'RM, 



(2) 



where C is the thermal capacity of the star. The expressions 
of C and L v for normal SQM are taken the same as those in 
1221 and 02511 . which were originally calculated by Iwamoto 
112 811 . In the presence of color superconductivity, the ther- 
mal capacity and neutrino emissions contributed by the paired 
quarks are significantly suppressed by some exponential fac- 
tors, which can be found in IU7I1 . 

The power of the DC heating is obviously determined by 
the transition rate of the nuclear matter at the bottom of the 
crust. Equation ([TJ further indicates that the transition rate 
depends on the variation rate of the spin frequency. Therefore, 
by denoting the heat release per dissolved neutron by q„, the 
power of the DC heating can be estimated by [22] 



Hoc 



-q n 



1 dM c dQ. 
nib dQ. dt ' 



where mj, and Q. are the mass of baryon and spin angular 
velocity of the star, respectively. In the conventional MDR 
model, the temporal derivation of the angular velocity can be 
calculated by dQ./dt - -Q/r m , where the magnetic braking 
timescale reads r m = 1.69 x 10 9 B^|(Q/ ^nGpY 2 s, fin is the 
magnetic field intensity in units of 10 12 G, and -\JnGp repre- 
sents the order of the Keplerian rotation limit with G being the 
gravitational constant and p the mean stellar density. How- 
ever, by considering the arising of r-mode instability, which 
can lead to the loss of stellar angular momentum via GR, 
the spin-down history of SSs needs to be reinvestigated more 
carefully. 

It is natural to consider that a great amount of the stellar ro- 
tational energy can be transferred into r-modes. Subsequently, 
the energy carried by the r-modes can be further dissipated by 
both bulk and shear viscosities. To be specific, the part of the 
r-mode energy dissipated by the bulk viscosity escapes from 
the star via neutrino emission, whereas the part corresponding 
to the shear viscosity can be deposited into internal energy, 
which could be of great importance for the stellar thermal evo- 
lution. Following [23] and [25], the power of such RM heating 
can be written as 



Hum — 



IE a 2 JMR 2 Cl 2 



(4) 



where E - ^a ■ JMR Q is the canonical energy of the r- 
modes with 7=1 .635 x 10~ 2 for n — 1 polytrope, a the ampli- 
tude of the r-modes, and M the stellar mass. For normal SQM, 
the shear viscous time scale is r sv = 5.41 x 10 9 a './ n' s, 
where a c ,o.i = ar c /0.1 is the strong coupling [8l I29I1 . In 
presence of color superconductivity, a fraction or all of the 
quarks are paired and thus the shear viscosity contributed 
by these paired quarks should be suppressed by a factor of 
exp(A/3kBT), where A is the pairing gap. In such a case, 
the shear viscosity due to electron-electron scattering be- 
comes more important for the r-mode damping. The cor- 
responding time scale can be expressed to t" = 2.95 x 
10 7 (// e //"q) 14 ^ 3 ?y s, where p. e and p. q are the chemical po- 
tentials of electrons and quarks, respectively. On the other 
hand, as pointed out by [8], an extra viscosity due to the rub- 
bing between the solid nuclear crust and the electron atmo- 
sphere of the quark core could lead to a more effective dissi- 
pation of the r-mode energy, which determines a time scale of 
T sr = 1.42 x 10 8 (v/l kHz)- I/2 r 9 s. 



III. SPIN DOWN OF SSS 

As mentioned above, the GR due to r-mode instability could 
play an important role in the spin down of a SS. Following the 
phenomenological model developed by Owen et al. 02311 and 
Ho and Lai 1 30], the spin down of the star affected by both 
MDR and GR brakings can be solved from 



dQ. 

Hi 



Q. 2Qa 2 Q. 



(5) 



(3) 



where Q = 0.094 is a constant determined by the stellar struc- 
ture. The viscous time scale here, r v , includes the effects of 
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FIG. 1 : Spin evolutions of NSSs (left panel) and CSSs (right panel). The solid and dotted lines correspond to the case of GR+MDR braking 
and of only MDR braking, respectively. Two different saturation values of k = 10~ 6 and 10~ 8 are respectively represented by the thick and thin 
lines. The typical magnetic fields of canonical pulsars of 10" G and 10 12 G are taken as labeled. 



both bulk and shear viscosities, which can be calculated by respectively, where k - a 2 at . As a rough estimation, the up- 



. For normal SQM, the bulk viscous time 
0.886(Q/ yftGpy 2 Tr 2 ™- 4 



m sjoo s ' where ot,,.i nn 



T v = (^bv+^v) 

scale reads Tbv = "•"'">--/ v"r» 

is the mass of strange quark in units of 100 MeV QSj, [29|]. In 
presence of color superconductivity, this time scale should 
also be prolonged by a factor of exp(2A/kBT). In order to 
solve Eq. ©, we must simultaneously determine the evolu- 
tion of the r-mode amplitude by 



da 
dt 



1 1 - a 2 Q 1 
- + =- 



(6) 



s for 



where the GR time scale reads r g = 3.26(Q/ ynG 
n — 1 polytrope. 

It is generally believed that r-mode instability would finally 
arrive at a saturation state due to some nonlinear effects. Lind- 
blom et al. J3ll |32[] suggested that the nonlinear saturation 
may be determined by dissipation of energy in the produc- 
tion of shock waves. Moreover, the decay of the amplitude 
of the order unity is due to leaking of energy into other fluid 
modes, leading to a differential rotation configuration 113311 . 
Afterwards, the coupling between r-modes and other modes 
was analyzed 1341 , 13511 . On the other hand, the role of differ- 
ential rotation in the evolution of r-modes was also studied 
thoroughly 11361 [3711 . All these works obtained a credible con- 
clusion that the maximal saturation amplitude of r-modes may 
not be larger than the small value of 10~ 3 . 

During the saturation state, Eqs. (0 and © would be re- 
placed by 



dQ. Q 1 2Q kQ 



and 



dt T m 1 — kQ t, 1- kQ 



da 



per limit of the r-mode energy can be connected to the stellar 
rotational energy by £ max =* KQE rot . 



In Fig. Q] we plot the spin-down behaviors of both NSSs 
and CSSs. By considering that only a small fraction of quarks 
in the 2SC phase are paired, the evolution behavior of a 2SS is 
qualitatively similar to a NSS (see [22]). So in this paper the 
case of 2SSs will not be discussed specifically. As previously 
found in 12511 . the spin velocity of SSs at the very first time 
could decrease exponentially due to the exponential increase 
of the r-modes. After the r-modes enter into the saturation 
state with a constant amplitude, the spin evolution dominated 



by the GR braking can be determined by Q oc -Q/t» 



-n 7 , 



which gives Q oc r I//6 . This analytical estimation could be 
pretty valid for CSSs, which is similar to the case of neutron 
stars Ir23ll24ll . but not so good for NSSs for t > 100 yrs due to 
a self-decrease of the saturation amplitude (Instead, Q oc r 1 ^ 14 
approximatively). After some time, the braking due to MDR 
would become comparable to and eventually exceed the GR 
braking, which leads to Q oc r 1 ^ 2 . In Fig. [1] two values of the 
saturation of k — 10~ 6 and 10~ 8 are considered for a compar- 
ison. Obviously, by pushing the saturation to smaller values, 
the influence of the r-modes becomes weaker. Then, the spin 
evolution returns to the traditional one due to only MDR and 
all issues discussed in this paper return to the case considered 
in Il22ll . So in the following calculations we only taken the 
value k - 10~ 6 into account. The other model parameters are 

(7) taken as follows. The initial values of the r-mode amplitude, 
angular velocity, and interior temperature are ao = 10" 10 , 
Qo = | yJTtGp, and Tq - 10 10 K, respectively. The param- 
eters for SSs are M - 1.4M , and R - 10 km, and for the 
SQM are a c = 0.2, A = 100 MeV, m s = 100 MeV, q„ ~ 20 

( 8 ) MeV, and Y e = 10~ 5 . 
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FIG. 2: The temporal dependences of the powers of the DC (solid) and RM (dashed) heatings for NSSs (left panel) and CSSs (right panel). 
The magnetic fields are the same as in Fig. Q] 
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FIG. 3: The evolutions of the surface temperature of NSSs (left panel) and CSSs (right panel) with both DC and RM heating effects. The 
magnetic fields are the same as in Fig. [T] The light grey bands present the case without heating. 



IV. RESULTS 



The temporal dependences of the heating powers are nu- 
merically calculated and presented in Fig. [2] As shown, the 
RM heating could paly a dominate role in the thermal evo- 
lution of SSs during their middle ages, i.e., from ~ 100 yrs 
to ~ 10 4 ~ 6 yrs. By contrast, the DC heating exceeds the RM 
heating for t < 100 yrs or t > 10 4 ~ 6 yrs. The stronger the 
dipole magnetic field, the weaker the RM heating. For NSSs 
with B = 10 12 G, the RM heating becomes always subordi- 
nate to the DC heating. On one hand, following the analytical 
spin-down history of CSSs, the temporal evolution of the DC 
heating can be approximated to Hoc K &■&■ K f 3 in the 
GR braking stage and Hoc K QQ K t~ 2 in the MDR braking 
stage, while for the RM heating H RM oc Q 2 /r sl oc r 5/12 T~ l 
and oc r 5 ^ 4 T~ l in the GR and MDR braking stages, respec- 
tively. On the other hand, for NSSs, we can get Hoc K f 8 ^ 7 
and following Hoc K t~ 2 , while Hrm oc £2 2 /t sv oc r lll T~ 5 ^ . 



The RM heating is eventually switched off due to the end of 
the r-mode instability. As a result, we plot the thermal evolu- 
tion curves of both NSSs and CSSs in Fig. [3] where the stan- 
dard cooling of SSs without heating is shown by the light grey 
band. As seen, the change in the thermal history of SSs aris- 
ing from heating effects could be remarkable, as previously 
claimed by [22]. 



For CSSs, as analyzed by Yu and Zheng [22], the main ther- 
mal evolution can be directly determined by the equilibrium 
equation L y - H, since both L v and C vanish due to the color 
superconductivity. Therefore, their thermal history during the 
equilibrium can be divided into RM and DC heating stages. In 
the first stage, T 2 oc r 5/12 r _1 (GR braking) and T 2 oc r 5/4 T _1 
(MDR braking) gives T oc t~ 5 ' 36 and T oc r 5 ^ 12 , respectively. 
The realistic one could fall between these two extremes. In 
the second stage, the DC heating and MDR braking gives 
T oc r 1 ^ 2 . Here, the most interesting thing could be the sharp 
transition between these two stages. As clearly shown in the 
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FIG. 4: The evolution of NSSs (left panel) and CSSs (right panel) for a ydt = 10 3 . The shaded areas are r-mode unstable regions. The labels 
for the open cycles represent the value of log(r/yrs). The magnetic fields are the same as in Fig. [TJ 



right panel of Fig. [3] an extremely steep decay of the stellar 
temperature appears in the thermal evolution curves of CSSs 
at the age of ~ 10 4 ~ 6 yrs. The range of the temperature dur- 
ing this transition is ~ 3 x 10 5 K— 3 x 10 6 K. Such a time and 
temperature region in the T — t panel is just consistent with 
the most observational data 1381 . On the other hand, such a 
steep decay indicates that the stellar temperature can be found 
to decrease, most rapidly, by about seven percent in ten years, 
as shown in the insert in the right panel of Fig. [3] Another 
interesting feature in the thermal evolution curves of CSSs is 
the big bump appearing during the first hundred years, where 
the equilibrium between the heating and cooling effects has 
not been built (.Hoc :s> L 7 ). Because of the appearance of 
the r-modes, the deconfinement transition of the stellar crust 
can take place much earlier and more quickly than the case 
with only MDR braking, which means that most binding en- 
ergy of the crust could be released within the first hundred 
years. Therefore, the bumps obtained here are nearly inde- 
pendent of the strength of the magnetic fields, whereas in 112211 
such a bump only exists in the high magnetic field case. It 
should be pointed out that, due to the thermal relaxation of 
the stellar crust which is not considered in our calculations, 
the variation of the interior temperature of SSs at early ages 
actually cannot be immediately reflected by the surface emis- 
sion. However, because the crust of SSs is only tiny, which 
is just the outer part of the crust of traditional neutron stars, 
the relaxation time scale of SSs could only be on the order of 
several years (e.g., see Figs. 2 and 3 in 1U6I0 . Therefore, an 
infant CSS at the age of several tens to a hundred years is still 
expected to have a high temperature of ~ 10 7 K (emitting at 
~ 1 keV). 

For NSSs, their thermal histories are not very different from 
the standard one, except for very late ages. In principle, the 
histories of NSSs could be divided into four stages, which 
are determined by Cf = H — L v , Cf = L v , Cf = L y , and 
CT - Hue (MDR braking), respectively. Of course, for a spe- 



cific NSS with a fixed magnetic field, there may only be two or 
three stages. In the first stage, the heating powers and the neu- 
trino luminosity could be comparable to each other, so an ana- 
lytical calculation is infeasible. For the following three stages, 
we can easily get T oc t I , T oc exp(-f), and T oc t > 2 , re- 
spectively, for C oc T and L v oc T 6 (QDU). Different from 
CSSs, it is not easy to produce an early temperature bump by 
NSSs because of their very intense neutrino emission during 
the first hundred years. 

Finally, for an integrated impression of the thermospin evo- 
lution of SSs, we plot the evolution curves of SSs in the T — v 
panel in Fig. |U where the r-mode instability windows are also 
presented. Some qualitative descriptions for such figures are 
provided in 12511 . Here we emphasize that (i) for NSSs, after 
their birth at the region with high v and T where the r-mode 
is stable, their thermospin evolution is basically determined 
by the boundary of the r-mode instability window and that (ii) 
for CSSs, they can stay at the bottom of the r-mode instability 
window for tens of thousands of years. 



V. SUMMARY AND DISCUSSIONS 

The SS hypothesis and the specific state of the SQM in SSs 
are some of the most fascinating mysteries in both physics 
and astronomy. In this paper, we investigated the thermal evo- 
lution of SSs by considering both the deconfinement transi- 
tion of the stellar crust and the r-mode instability of the fluid 
core. The results show the following: (i) The cooling of NSSs 
with heating effects is only somewhat different from the well- 
known standard one. Therefore, there seems still no signifi- 
cant discrepancy between the NSS model and the temperature 
observations of pulsars, (ii) For CSSs, there are two interest- 
ing features found in their thermal histories. One is the big 
temperature bump during the first hundred years and the other 
is the steep decay of the temperature at the ages of ~ 10 4 -10 6 



yrs. 

The serious problem of the CSS model is that, to our knowl- 
edge, neither the early bump nor the steep decay (~7% in 
ten years) has been reported from observations, though they 
should be easily detected if they indeed exist. Let us recall 
the constraint on the CSS model given by Madsen |j§], who 
claimed that the CSS model is not permitted by the most rapid 
pulsars because these pulsars locate deep at the r-mode in- 
stability window of CSSs. Therefore, we conclude that both 
the rotational and thermal tests on the SS models, which are 
independent of each other, suggest that the CSS model is dis- 
favored by both the rotation and temperature observations of 
compact stars. As an alternative consideration discussed in 
18[], the following two possibilities are worthy of attention: (i) 
The CFL phase is not reached at densities relevant in pulsars, 
but the SQM could still be stable because the NSS model is 
not ruled out. (ii) The CFL phase appears in a hybrid star, in 



which only the stellar core consists of quark matter. Similar to 
the traditional neutron stars, hybrid stars could pass the tests 
by both the spin and temperature observations of pulsars. 
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